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Abstract 



We present a simple, exact and self-consistent cosmology with a phenomenological model 
of quantum creation of radiation due to decay of the scalar field. The decay drives a non- 
isentropic inflationary epoch, which exits smoothly to the radiation era, without reheating. 
The initial vacuum for radiation is a regular Minkowski vacuum. The created radiation obeys 
standard thermodynamic laws, and the total entropy produced is consistent with the accepted 
value. We analyze the difference between the present model and the model with the decaying 
cosmological constant considered in Q. 
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1 Introduction 



The aim of this paper is to extend the scenario of the evolution of the universe with smooth 
exit from inflation, and particle production at the expense of the decaying cosmological 
constant, developed in paper IQ. In that paper thermodynamics and Einstein's equations led 
to an equation in which the Hubble rate H is determined by the particle number N. The 
model is completed by specifying the particle creation rate F = N/N, which led to a second- 
order evolution equation for H. The evolution equation for H then has a remarkably simple 
exact solution, in which a non-adiabatic inflationary era exits smoothly to the radiation era, 
without a reheating transition. For this solution, there were given exact expressions for 
the cosmic scale factor, energy density of radiation and vacuum, temperature, entropy and 
super-horizon scalar perturbations. 

In this paper we would like to generalize the abovementioned results for the case of a 
scalar field ip interacting with radiation via the gravitational field, leading to cosmological 
evolution with smooth exit from inflation. Our particular task is to determine whether the 
theory formulated in terms of the scalar field 93 leads to any new results in comparison with 
the previous case of a decaying cosmological constant. 

The background of this paper is constituted by two ideas in physical cosmology. One 
of them is connected with the longstanding attempt to explain all matter in the universe as 
produced by quantum creation from vacuum. This has been studied via quantum field theory 
in curved spacetime (see for example [§, Most cosmological models exhibit a singularity 
which presents difficulties for interpreting quantum effects, because all macroscopic param- 
eters of created particles are infinite there. This leads to the problem of the initial vacuum 
(see im ) . One attempt to overcome these problems is via incorporating the effect of particle 
creation into Einstein's field equations. For example, in the papers of the Brussels group |^, 
the quantum effect of particle creation is considered in the context of the thermodynamics of 
open systems, where it is interpreted as an additional negative pressure, which emerges from 
a re- interpretation of the energy-momentum tensor. This effect is irreversible in the sense 
that spacetime can produce matter, leading to growth of entropy, while the reverse process 
is thermodynamically forbidden. 

The main difference with our present paper is that in Q the law of massive particle 
creation, i.e. the mechanism of energy fiow from the gravitational field to matter, leads to 
a non-zero number of particles at the beginning of expansion, described as a fluctuation of 
the regular vacuum. These results were recently generalized in a covariant form in Our 
approach differs from that of [§, ^ in that we do not modify the field equations. Instead, we 
associate the source of created particles as a decaying vacuum of the infiaton field (p. 

A number of decaying vacuum models has appeared in the literature (see 0, ^, ^ and 
references cited there). A review of the different phenomenological models of evolution with 
variable cosmological term can be found in paper ]lO| . 

Infiationary models with fixed cosmological constant and cold dark matter have been 
successful in accounting for the microwave background and large-scale structure observations, 
while also solving the age problem (see [11 1). However, these models are challenged by the 
reduced upper limits on A arising from the Supernova Cosmology Project, and also by the 
long-standing problem of reconciling the very large early-universe vacuum energy density 
with the very low late- universe limits [^. 

One resolution of these problems is a decaying cosmological constant A which is treated 
as a dynamical parameter. This approach was typical for the quantum field theorists for 
many years (see for example [|l2|). Anything which contributes to the energy density 
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of the vacuum behaves like a cosmological term via = SirGpy. Many potential sources 
of fluctuating vacuum energy have been identified (including scalar fields which were 
to give rise to a negative energy density which grows with time, tending to cancel out any 
pre-existing positive cosmological term and drive the net value of A toward zero. Processes 
of this kind are among the most promising ways to resolve the longstanding cosmological 
'constant' problem (see |14] for review). It is worth mentioning the recent paper of Parker 
|15| indicating an attempt to revive the idea of the cosmological constant as a purely quantum 
effect associated with the renormalization of the general relativistic action. 

In ad hoc prescriptions, the functional form of A(t) or A(a) or A{H) (where a is the scale 
factor and H is the hubble rate) is effectively assumed a-priori (see the review ||lO|| where all 
known forms for A are listed). Typically, the solutions arising from ad hoc prescriptions for 
A are rather complicated, and moreover, it is often difficult to provide a consistent simple 
interpretation of the features of particle creation, entropy and thermodynamics. 

In common with Q and extending Q, we attempt to provide some clear and consistent 
physical motivation for the particular form of vacuum decay for the field ip. In contrast to 
many other models, we propose a simple, exact and thermodynamically consistent cosmo- 
logical history. The latter originates from a regular initial vacuum. Together with naturally 
defined asymptotic conditions for the number of created particles and the form of the poten- 
tial for the field <p this leads to a simple expansion law and thermodynamic properties, and 
to a definite estimate for the total entropy in the universe. 

Non-adiabatic inflationary models differ from the standard models (see for example ||l^ ) , 
in that: (a) radiation is created continuously during inflation, rather than during reheating; 
(b) the continuous vacuum decay itself initiates a smooth exit from inflation to the radiation 
era; (c) entropy and heat production take place continuously, without the need for reheating. 
In the standard approach, the scalar field drives adiabatic (i.e., isentropic) inflation, followed 
by a non-equilibrium reheating era when the field decays into radiation and inflation is brought 
to an end. The potential of the field is then the key ingredient. We argue in this paper that 
the potential (treated usually as a self-interacting potential) actually consists of two parts. 
One of them is the interaction of the field with radiation via the gravitational field. This 
term is expressed in purely geometrical terms and in the case when the field is zero plays the 
role of the decaying cosmological constant as in paper [Q. The second part of the potential, 
corresponding to self-interaction of the field, is postulated in a simple form consistent with 
the standard requirements of inflationary cosmology. One should mention that our ansatz 
of decomposing the potential into two parts is in some sense similar to the postulate of the 
interaction term in the form T ip^ which was used e.g. in the papers ||8|, In these 

papers the strategy was to modify the second order Klein-Gordon equation. In our case the 
decomposition of the potential allows one to formulate a differential equation of the first order 
for the field ip. 

Since the exit from inflation to the radiation era is smooth, we avoid the problem of 
matching at the transition. A similar smooth evolution has been used in |19, but in the 

context of adiabatic inflation, and without a consistent physical foundation. In analogy with 
Q we show that the ad hoc form of H{a) given in [19| follows from simple physical conditions. 
In 1^], a kinematic analysis is given for various non-adiabatic inflationary evolutions with 
smooth exit, but these evolutions are outside the scope of our model. 

The choice of a as dynamical variable and the very simple form of H{a) that meets the 
physical conditions, lead to elegant expressions for all parameters describing the radiation 
and decaying vacuum, and also to a physically transparent interpretation of these results, 
including the estimate of entropy. 
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We use units with SvrG, c and A;b equal 1. 



2 Model 

Consider a spatially flat FRW universe 

ds"^ = -dt^ + a'^{t) dx^ + dy^ + dz^ , (1) 

containing radiation with energy density pr{t) and pressure pr = Pr/'i- 
The energy momentum tensor of these components correspondingly is 



We also consider quantum matter described by the scalar field ip with energy density and 
pressure given in a conventional form 

where ^ is a potential. In a special case of quantum matter with the equation of state 
Pip = — Pip , we have ip'^ = 0, and the potential can be treated as a cosmological constant, 
dependent on t. 



3 Conservation Law 

The conservation equations V^{T^j^ + T^) = in the general case for 7^ reduce to 

Pr + 4:Hpr = - p^ - 3H if^ (2) 

where H = d/a is the Hubble rate. The equation (^) shows how energy is transferred from 
the quantum field to radiation. This energy transfer can be understood as creation of quanta 
of radiation at the expense of the decaying field Employing the extended form of the first 
law of thermodynamics suggested in [Q] and applied to radiation only: 

d( Pr V) + PrdV - —d{nrV) = , (3) 
rir 

where hr = pr +Pr = 4/3 pr , one can link the evolution of pr with the change of the total 
number of photons Nr = UrV, where V is the comoving volume of the observable universe. 
Since (^) is equivalent to 

4 Nr 

Pr +4H Pr = - Pr (4) 

a comparison with (^) gives for the particle creation rate 

Nr _ Pip + 3H ip^ 



Nr 4/3 Pr 



(5) 



The creation of radiation implies that Nr/Nr > 0, which means in turn that subject to 
Pr > 0, one must have pip + 3H ip^ < 0, i.e. the energy of the scalar field ip decreases in 
time, being converted thereby into quanta of radiation. 
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In the case ip^ =0, p^p =V = h.so that the equation @ becomes 



(this case is considered in i.e. the production of classical photons is a result of decay of 
the vacuum with the energy density A. 

If if^ =0 and V = const (or A = const) the initial vacuum for photons (where A^^ = 0) 
is stable leading to no particle production. The total number of particles is conserved and 
equals zero. 

By switching on the source of radiation p^p 7^ in the right-hand side of the equations 
(|2|) and (^) we effectively switch on the coupling of radiation with gravitational field, leading 
to creation of photons from vacuum. 



4 Field equations 

The field equations 



R. 



are 



3H^ 



rpR _|_ rpQ 



2H + 3H^ 



Pr + Pip 
1 



9^' +V, 



(7) 
(8) 



and if both are satisfied then the energy equation (^) follows identically. 
Combining (^) and (^ one finds 

Pr =-^[2H+ ip^ ], 



(9) 
(10) 



where R = QH + 12H^ is the Ricci curvature. It is easy to see that the equation dl 



equivalent to the equation linking the potential V and (p"^ : 



IS 



V 



R 



4^ 2 4^ 4 



(11) 



We can use the equation (|5|) in order to express Nr in terms of pr , and then in terms of the 
combination 2H + ip"^ : 



C- 



_iV4/3 



N = C 



-{2H+ p^ )a 



l3/4 



(12) 



The formulas (^, 10, |T2|) give the expressions for the physical quantities describing matter 
in the universe in terms of H and <p^ , which have to be determined from two additional 
conditions. First we need the equation for (p and then we need another constraint in order 
to define H. 

The standard procedure to find 99, and hence p^ , is to obtain the analogue of the Klein - 
Gordon equation for the field v?. Indeed, substituting pr @ and (|lO|), into the equation 
(|5|) one can get the equation 



ip ip + V + 3H ip^ =[2H + ip'^ ]r. 



(13) 
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The difference between tliis equation and the conventional form of the Klein-Gordon equation 
(see for example 18 1) is that we can not cancel ip and proceed to an equation of second 
order for (/? because of the geometrical term 2HT at the right hand side of this equation. 
Using the expression for V (O) the last equation can be rearranged in the following form 



-(p p +[3H 



r] p 



•2 



4 



(14) 



From here one can easily realize that the standard Klein-Gordon equation for ip will be 
achievable only in the case when the right hand side of ( |l^ ) equals zero. 

It is interesting to see that in the case 99^ = the right hand side of the last equation 
will be zero anyway, giving the connection formula between the geometrical parameters of 
the model and the rate of particle creation: 



2HT{a) - V = 
4 



(15) 



This equation has already been obtained in our previous model |j^] (equation (9)). It was 
solved there with respect to iif as a function of a under some physical assumptions about 
the rate of particle creation F. This led to a scenario with a smooth exit from inflation to 
radiation dominated era (see details in [Q]). In the case of the paper Q it was sufficient to 
impose one physical condition for F in order to determine H{a). 

It is interesting however in the present situation that one can get a nontrivial solution for 
99 7^ if the condition (15) holds. In this case we have 



dt 



+ 



AH 



R 
6H 



0. 



(16) 



Following P, |, we use a as a dynamic variable instead of t, and consider the Hubble 
rate as H = H{a) (in this case we can not consider a = constant as a limiting case for a flat 
universe). Rewriting the last equation in terms of a, and using the prime as a notation for 
the derivative with respect to o {d/dt = aHd/da = aH{)') one can rewrite the equation ( [T^ ) 
in the form 



da 



(17) 



This equation integrates to 



A{-aHH') = A{-H) 
Substituting this result into the formula for pr one gets 



Pr 



A'] i-H) 



which implies that A"^ < 3/2. This formula together with the condition ([T5| ) reproduces the 
result of the paper Q . 

Coming back to the general case of the equation (|l^, 

Z + P{H,T)Z = Q{H,T) (18) 

where Z = (p"^ , P = A H — A/SF, Q = A/3{2HT — R/A), one can simplify it using the following 
substitution 

Z = Zh — 2H, 
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which gives the equation for (compare with the equation (18), when Q = 0): 



Z/, + p(i/,r)z^ = o 

This equation can be easily solved in terms of N: 

Zh{a) = CN^/^{a)a-^, so that Z = ^'^ = -2H + CN^/^{a)a-^ 

(where C is an arbitrary constant) which is entirely equivalent to (|l2|). In other words in 
order to determine if"^ one needs to know N and H, and this shows that we have exhausted 
all information about our system from the conservation laws (|2|), dH) and the field equation 
(|^. One should develop some further arguments in order to determine 0^ and H. 



5 The potential 

The expression for the potential V which we obtained above (^) is very interesting because 
it consists of two parts, one proportional to the kinetic energy of the field (p"^ , another 
explicitly a function of time or a, which describes the change of energy in the system field- 
radiation due to expansion of the universe. This makes it possible to argue that the potential 
can be interpreted in the following way: 

V{ip,a) = Uiip) + Vint{a), (19) 

where 

U{ip) = \^\ VUa) = ^ (20) 

The potential U describes the self-interaction potential of the field (p, whereas Vint corre- 
sponds to interaction of the field (p with radiation via the gravitational field. In the case of 
no gravitational field, i.e. H = 0, V = 0. 

Now we are in a position to solve equation ( [20| ) with the initial condition f7 — > as a — > 
(which follows from (^) as an initial condition for ip'^ (a = 0) = 0). The solution depends, 
obviously, on the form of U. The contribution of p"^ in the formula for pr can be 
treated in this case as an additional external source of radiation. Solving the equation ( [20| ) 
does not give, however, any information about H{a), so that we still need some additional 
physical arguments in order to fix H{a). 

The choice of U is not a trivial one in our case. Since we would like to develop a scenario 
of the evolution of the universe with infiation it is natural to assume that the potential can 
be chosen in a conventional form 

C/ = t/o-imV + ^^^ 

which implies the so called 'rolling down' of the field from its unstable value 93 = with 
U = Uq to stable asymptotic value ipo which is defined from the equation U{ipo) = 0. This 
potential, however, must explicitly satisfy the 'initial' condition U — > as a ^ 0. Assuming 
that the rolling down of the field if starts at a = we find ourselves in a difficult situation 
because the potential C/ as a function of f must have a nonzero value for the initial value 
(p = whereas the same potential must be equal zero at a = as a function of a. 

One possible solution of this problem is to take into account that the potential U after 



solving the equation (20) is a function of a. Since we know a priori what properties of the 



model to expect we can postulate U in the form 

U = a'^H^\{p^ - ^If (21) 
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where A and (fo are constants. The theory with this potential will also benefit from the fact 
that the dynamical equation (|2^) can be solved in terms of the variable a without knowledge 
ofH{a). 

This potential satisfies naturally the condition U ^ as a ^ 0. From one can 
intuitively conclude that the field evolves to i^po which is the turning point of the potential. 
The dimensionless constant (fo is a free parameter and corresponds to the value of the field at 
stable vacuum with zero energy. It is clear that the potential (21) has a maximum at (/? = 0, 
so that one can argue that the value of Uq = Xip^ determines a temporal scale (in terms of 
a) of decay of the unstable vacuum into particles as it is in a standard inflationary scenario. 
We show below that the parameter [ly/Vo) ^ can be interpreted geometrically as the value 
of the scale factor at the point of exit from inflation. 

The equation (]20|) now has the form 



2\2 



and its solution is 



1 + C exp ib4\/A(^oQ^ 
l-C exp ±4\/A(^, 



From the asymptotic condition (/9 — > (/?o as a — > cxd we choose the sign (— ) in the exponential 
functions; then from the initial condition 99(0 = 0) =0 one finds C = — 1 so that finally 



93 = 



1 — exp(— 4\/A99, 



1 + exp(— 4\/Av'o'^) 



tanh[2 \/A 



ifQ a\ 



(22) 



This solution describes a standard rolling down of the field if along the potential curve 
U from the unstable value (/? = at a = (where 99'^ / 0) to the stable state with ip = i^q 
and zero potential. From (^) one can find 



' a ' 


2 


'2 a' 




cosh~^ 








. (/Jo a* . 



(23) 



where we introduced a parameter a* = {(pgy/X)"^. The field f evolves effectively to its stable 
value ifo during the interval (a*(/?o)/2. This implies that the contribution of ip"^ to pr @ 
and N (|l2|) is exponentially small for a > a*(/3o- 



6 The Hubble parameter H(a) 

Now we are in a position to make a prediction of the asymptotic form of the Hubble parameter 
Hia). 

Making a natural physical assumption that N{a) A^oo < 00 as a ^ cx) (in analogy 
with Q), one can obtain from (|l2|) the asymptotic for H{a) as a ^ a^ipQ neglecting the 
contribution from (p"^ in the formula for N (|T2|): 



H{a) ~ a ^ as a ^ 00. (24) 
In the region a <^ {a^:p>)/2 one can use an approximate expression for ip'^ 



a 
a* 



2 
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which, being substituted into the equation (^) leads to the differential equation for H^{a): 



4a 



[H\a)]' + ^H\a) = -C 



(25) 



In the limiting case N{a) 
equation 



as a — > one can obtain a formal nontrivial solution of this 

2 



H^{a) ~ D exp(-2 



which gives an accurate asymptotic for H as a — > 0: 

2\ 



H'^(a) ^ D ll 



a 
a* 



as 



0. 



(26) 



Comparing ( p4| ) and (Pq) one can argue that the simplest smooth function H{a) satisfying 
both these asymptotic conditions is 



H{a) 



C 



ai + 



where C, B are constants. This describes a universe with a smooth transition from inflation 
to radiation. The universe is initially de Sitter-like (since H ~ constant for small a), and 
becomes radiationlike (since ~ 1/a^ for large a) . The parameter a* can have a geometrical 
interpretation now as a value of the scale factor at exit from inflation, defined in general 
geometrically from the condition d(tc) = 0, or equivalently 



-aeH'^, i.e. 



a{te). 



Thus the form of the Hubble rate is 

H(a) 



2H, 



" a2 + o2 



(27) 



This form was presented in as an ad hoc prescription to achieve smooth exit from inflation 
to radiation, but without a physical basis such as that given here. In the paper j^] we obtained 
the same form of the Hubble rate in the model with a smooth exit from inflation as an exact 
solution of a differential equation subject to a hypothesis on the rate of particle creation in 
the model. The novelty of this paper is that it is possible to link the parameter of exit from 
inflation Oe with the parameters of the self-interacting potential of the scalar field f. 

ae = a, = {^lVX)-\ (28) 

The expression for the cosmic proper time follows on integrating Eq. ([27| ) (see also |j^): 



te + 



AH, 



ln(^ 



+ (-| -1 
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7 Thermodynamics of radiation 



On substituting now the expression for (jp' ( p3[ ) with a^, 
density of radiation (^), one obtains 



where 



Pr {a) = Pria) 



p°(a) = 12H^, 



ttg into the formula for the energy 



-V'o , (29) 



is the energy density of radiation in the case of ip = (compare with 



and 



F 



1 + 



cosh 



2 a 



(30) 



(31) 



From the physical condition pr > 0, and hence from F < 1, one can find the range of values 
for ifo: 

0<ipo< 2V2, (32) 

which in conjunction with (^) determines the range of possible values of A (subject to our 
knowledge of Oc). It is clear that F{a = 0) = 1. It follows from (p2[) that for all possible 



values of (po from (|32D the field ip decays on scales a < Og, so that it is impossible to extend 
the decay of for a > Cg. 

For the energy density of the field ip we have 



(a) = A(a) 



1+ - F -,^0 



where 



A(a) 



12Hf 



(33) 

(34) 
It follows 



is the Pip when = i.e. the decaying cosmological constant of the paper [ 
that p^{0) = 12H^. Note that (||) implies H{0) = 2He. 

The function F represents a rapid exponential decay so that one can effectively treat it 
as zero for a > Cg. This gives the expected asymptotics for pr and p^p as a — > 00: 



p{a) 



Pf (a) 



so that pp: rapidly becomes negligible in comparison with p^ ■ 

Comparing the functions p^, Pr Pip in terms of the dimensionless variable x = aja^ 
one can conclude that the contribution of 93^ in the expressions for pr and p^ , (which 



appears in the formulas ( ^91 ) and (33) through the function F), is noticeable only for a < a^. 
The order of magnitude for p^ and pj?, and for p^p and A correspondingly is the same, 
leading only to some shift of the point of maximum for towards exit from inflation. Note, 
that reaches a maximum at Om = ac/\/2, with 



P?n^P°(«m) = f ^e, A(a^) = 2/,^ 
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Note also from the graph that pr and A are equal at exit up to terms of 0(10 

Pr (ae) « Pr "(Oe) = = A^a^) « p^ (Ce) , 

while p^ <C Pr for a ^ a^, i.e. during the radiation-dominated era. 

The formulas (|29| ) and (^) reflect the creation of radiation due to vacuum decay. The 
initial value pr (0) = confirms that the field corresponding to radiation is initially in a 
regular vacuum state. 

Substituting the equation ( |27| ) into equation (|l2|) we get the exact form for the particle 
number 



N{a) = iV°(a) 

where 

iVO(a) = N, 



l-F{-,ipo 



l3/4 

(35) 



is the particle number in the case 99 = 0, and N^^o is usually taken to be about 10^*^ (see [|l|). 
It is clear that the presence of F in the formula ( |35|) does not affect the asymptotic behavior 
of as a ^ 00, dying away exponentially on scales of the order of Cg. It means that after 
exit from inflation the behavior of radiation and the decay of the energy density for (/? will 
be practically the same as they were in the case of the cosmological constant decay Q. 

Since A'^(O) = n(0) = 0, the initial state of the field has no particles, i.e. it is a regular 
vacuum. The number density is 



nia) 



X 3 / 2 \ 9/4 
cosh2 25/4ne(^) H-Ff/'. (36) 



In analogy with [T| it seems reasonable to use the black-body relation for the radiation 
throughout the expansion, and to define the temperature by 



1 p{a) i?223/4 



ap \ / a 



2 



3/4 



T{a) = -^= (^) I J (l-F)^/\ (37) 

^' 3n{a) ncCosh2 V a / I + aV ^ ^ ' ^ ^ 



where we have used equations ( p9D and (36). At the initial radiation vacuum, it is clear that 
T(0) = 0. During the radiation era, i.e. for a ^ Oe 



T ~ a 



-1 



in agreement with the standard result for free radiation in an expanding universe. 
The formulas for p and n can be presented in the thermodynamic form 



p=24 I T', n = 8 I T'. (38) 



Combining now the Gibbs equation 

TdS = d(pV) + pdV , 



with equation (|3|), and using the definition ( |37D of T, we obtain the entropy of radiation in 
the observable universe as 

S{a) = m{a) , (39) 
and leads to a value of the same order of magnitude as our result. 
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8 Conclusion 



One can conclude that the presence of the scalar field in this model does not change 
considerably the physical results which have been obtained in the paper Q with the decaying 
cosmological constant. The small difference with the paper ||l| can be observed only for 
a < Og. One can not extend the decay of the field ip beyond a = Oe in this model because 
of the condition on the range of values of (/Jq (32). This means that the thermodynamic 
analysis for o > Ce is similar to the paper |^] and gives the same prediction for the total 
entropy produced in the universe. 

Generalizing the results of this paper one can claim that models with decaying cosmolog- 
ical constant A corresponding to a special case of the equation of state p^p = — p^p ^ describe 
adequately the smooth transition from inflation to radiation and give a reasonable prediction 
for the entropy of matter in the universe. 
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